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Abstract 

The essentials of a new method in solving very large classes of nonlin- 
ear systems of PDEs, possibly associated with initial and/or boundary 
value problems, are presented. The PDEs can be defined by continu- 
ous, not necessarily smooth expressions, and the solutions obtained 
cab be assimilated with usual measurable functions, or even with 
Hausdorff continuous ones. The respective result sets aside completely, 
and with a large nonlinear margin, the celebrated 1957 impossibility 
of Hans Lewy regarding the nonexistence of solution in distributions 
of large classes of linear smooth coefficient PDEs. 



1. The Class of Nonlinear Systems of PDEs Solved 



The nonlinear systems of PDEs considered in this paper are composed 
of equations of the general form 

(1.1) F(x,U(x), . . . ,DPU(x), . . . ) = f{x), xett C R n 

were the domains Q can be any open, not necessarily bounded subsets 
of M n , while p G N™, \p\ < m, with the orders meffof the PDEs 
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arbitrary given. 

The functions F which define the left hand terms are only assumed to 
be jointly continuous in all of their arguments. The right hand terms 
/ are also required to be continuous only. 

However, with minimal modifications of the method, both F and / 
can have certain discontinuities as well, [3] . 

As it turns out, regardless of the above generality of the nonlinear 
systems of PDEs considered, and of possibly associated initial and/or 
boundary value problems, one can always find for them solutions U 
defined on the whole of the respective domains Q. These solutions 
U have the blanket, type independent, or universal regularity property 
that they can be assimilated with usual measurable, or even Hausdorff 
continuous functions, [1,4-6]. 

Thus when solving systems of nonlinear PDEs of the generality of 
those in (1.1), one can dispense with the various customary spaces of 
distributions, hyperfunctions, generalized functions, Sobolev spaces, 
and so on. Instead one can stay within the realms of usual functions. 
However, functional analytic methods can possibly be used in order 
to obtain further regularity or other desirable properties of such solu- 
tions. 



2. A Natural Smooth Framework for Nonlinear 
Systems of PDEs 

Let us now associate with a nonlinear PDE in (1.1) the corresponding 
nonlinear partial differential operator defined by the left hand side, 
namely 

(2.1) T(x, D)U(x) = F(x,U(x), . . . ,D%U(x), ...), xeVl 

Two facts about the nonlinear PDEs in (1.1) and the corresponding 
nonlinear partial differential operators T(x,D) in (2.1) are important 
and immediate : 

• The operators T(x, D) can naturally be seen as acting in the 
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classical or smooth context, namely 



(2.2) T(x,D) : C m {VL) 3 U 



T(x,D)U G C°(Q) 



while, unfortunately on the other hand : 

• The mappings in this natural classical context (2.2) are typically 
not surjective even in the case of linear T(x,D), and they are 
even less so in the general nonlinear case of (1.1). 

In other words, linear or nonlinear PDEs in (1.1) typically cannot be 
expected to have classical solutions U G C m (Q), for arbitrary contin- 
uous right hand terms / G C°(Q), as illustrated by a variety of well 
known examples, some of them rather simple ones, see [3, chap. 6]. 
Furthermore, it can often happen that nonclassical solutions do have 
a major applicative interest, thus they have to be sought out beyond 
the confines of the classical framework in (2.2). 

This is, therefore, how we are led to the necessity to consider gen- 
eralized solutions U for PDEs like those in (1.1), that is, solutions 
U C m (Q), which therefore are no longer classical. This means that 
the natural classical mappings (2.2) must in certain suitable ways be 
extended to commutative diagrams 

T(x,D) 

C m {Q) - C°(fi) 



(2.3) 



C 



c 



X 



Y 



T 



with the generalized solutions now being found as 



(2.4) UeX\C m (tt) 



instead of the classical ones U G C m (Q) which may easily fail to exist. 
A further important point is that one expects to reestablish certain 
kind of surjectivity type properties typically missing in (2.2), at least 
such as for instance 
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(2.5) C°{Q) C T(X) 



As it turns out, when constructing extensions of (2.2) given by com- 
mutative diagrams (2.3), we shall be interested in the following some- 
what larger spaces of piecewise smooth functions. For any integer 
< I < oo, we define 



3 T C O closed, nowhere dense 
u G C\tt\T) 



(2.6) C l Jfl) = | u : n -> 1 
and as an immediate strengthening of (2.2), we obviously obtain 

(2.7) T(x,D)CZ(Q) C C° d (Q) 

The solution of the nonlinear PDEs in (1.1) through the order com- 
pletion method will come from the construction of specific instances 
of the commutative diagrams (2.3), given by 



T(x,D) 



(2.8) M%(n) 



T 



bijective 



where the operation ( ) # means the order completion, [3], of the re- 
spective spaces, as well as the extension to such order completions of 
the respective mappings, see (5.7) below. It follows that in terms of 
(2.3), we have 

X = M^{Qf, Y = M°{Qf, T = T* 

thus we shall obtain for the nonlinear PDEs in (1.1) generalized solu- 
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tions 



(2.9) U e M^{Q)# 

Furthermore, instead of the surjectivity condition (2.5), we shall at 
least have the following stronger one 

(2.10) Ck(n) C T#(M^(Q)#) 

So far about the main ideas related to the existence of solutions of gen- 
eral nonlinear PDEs of the form (1.1). Further details can be found 
in [3,1,4-6]. 

As for the regularity of such solutions, we recall that, as shown in [1], 
one has the inclusions 

(2.11) M°(Q)* C Mes{Q) 

where Mes (Q) denotes the set of Lebesgue measurable functions on Q. 
In this way, in view of (2.8) and (2.9), one can assimilate the general- 
ized solutions U of the nonlinear PDEs in (1.1) with usual measurable 
functions in Mes (Q). 

Recently, however, based on results in [1,4-6], it was shown that in- 
stead of (2.11), one has the much stronger regularity property 

(2.12) M°{tif C M(fi) 

where EI (Q) denotes the set of Hausdorff continuous functions on Q. 
Consequently, now one can significantly improve on the earlier regu- 
larity result, as one can assimilate the generalized solutions U of the 
nonlinear PDEs in (1.1) with usual functions in HI(fi). 

Regarding systems of nonlinear PDEs such as in (1.1), with possibly 
associated initial and/or boundary value problems, it was shown in [3, 
chap. 8] the way they can be dealt with the above order completion 
method. 

In this respect, a surprising advantage of the order completion method 
is the ease, when compared with the usual functional analytic ap- 
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proaches, in dealing with initial and/or boundary value problems. 

3. Solution Method by Order Completion 

As first introduced and developed in [3], the mentioned existence and 
blanket, type independent, or universal regularity property of solutions 
for nonlinear systems of PDEs of the above generality in (1.1), together 
with possibly associated initial and/or boundary value problems, is 
based on a rather simple, however, quite general order completion of 
suitable spaces of piecewise smooth functions. 

Here one can mention that, in fact, this order completion based solu- 
tion method reaches far beyond the solution of PDEs, and in fact it 
can be applied to the solution of general equations 

(3.1) T{A) = F 
where 

(3.2) T:X — > Y 

is any mapping, X is any nonvoid set, while (Y, <) is a partially or- 
dered set, or in short, poset, while F G Y is given, and A G X is the 
sought after solution. 

In general, for a given F G Y , there may not exist any solution A G X 
for the equation (3.2), unless the setup (3.1), (3.2) is extended in suit- 
able ways. 

The usual such extensions assume suitable topologies on X and Y , 
and certain continuity properties for the mapping T in (3.2). 
As shown in [3], and also seen in the sequel, for the same purpose 
of solving the equations (3.1) in a suitably extended setup, one can 
successfully use the order completion of the spaces X and Y . 

A major advantages of the order completion approach is that such a 
method 

• does no longer differentiate between linear and nonlinear opera- 
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tors T in (3.2), in case the spaces X and Y may happen to have 
a linear vector space structure, [3]. 

Two further convenient features of the order completion method in 
solving general equations (3.1) are the following : 

• one obtains necessary and sufficient conditions for the existence 
of solutions, 

• one obtains explicit expressions of the solutions, whenever they 
exist. 

The order completion method for solving general equations (3.1) pre- 
sented here is but one of the possible such approaches. This method 
is based on so called "pull-back" partial orders on the domains of the 
mappings T in (3.1), (3.2). However, the order completion method 
for solving general equations (3.1) is not limited to "pull-back" par- 
tial orders. Further possible developments in this regard of the order 
completion method, and no longer based on "pull-back" partial orders 
will be indicated. 

Here it should be mentioned that usual functional analytic methods 
for solving linear or nonlinear PDEs often make use of "pull-back" 
topologies on the domains of the respective partial differential opera- 
tors as illustrated next in section 4, see for further details in [3, chap. 
12]. 

4. Use of "Pull-Back" in Functional Analytic Solution 
Methods of PDEs 

We present here one of the classical examples where a "pull-back" 
topology on the domain of a partial differential operator is used in 
order to solve the corresponding PDE. Let us consider on a bounded 
Euclidean domain Q, which has a smooth boundary d Q, the follow- 
ing familiar linear boundary value problem, usually called the Poisson 
Problem 
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A U(x) = f(x), xett 

(4.1) 

U = on dQ 

As is well known, for every given / G C°°(Q), where Q denotes the 
closure of Q, this problem has a unique solution U in the space 

(4.2) X = { ^r(0) I v = on d Q } 
It follows that the mapping 

(4.3) X 3v i — ► || Av || L2( „ } 

defines a norm on the vector space X. Now let 

(4.4) Y = c°°(n) 

be endowed with the topology induced by L 2 (f2). Then in view of 
(4.1) - (4.4), the mapping 

(4.5) A : X -> Y 

is a uniform continuous linear bijection. Therefore, it can be extended 
in a unique manner to an isomorphism of Banach spaces 

(4.6) A : X -> Y = L 2 {Q) 

In this way one has the classical existence and uniqueness result 
V /Gl 2 (fl) : 

(4.7) 3 \ U eX : 

AU = f 

The power and simplicity - based on linearity and topological comple- 
tion of uniform spaces - of the above classical existence and unique- 
ness result is obvious. This power is illustrated by the fact that the 
set Y = L 2 (Vt) in which the right hand terms / in (4.1) can now be 
chosen is much larger than the original Y = C°°(Q). Furthermore, 
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the existence and uniqueness result in (4.7) does not need the a priori 
knowledge of the structure of the elements U G X, that is, of the 
respective generalized solutions. This structure which gives the regu- 
larity properties of such solutions can be obtained by a further detailed 
study of the respective differential operators defining the PDEs under 
consideration, in this case, the Laplacean A. And in the above specific 
instance we obtain 

(4.8) X = H 2 {Q)nH^{Q) 

As seen above, typically for the functional analytic methods, the gen- 
eralized solutions are obtained in topological completions of vector 
spaces of usual functions. And such completions, like for instance the 
various Sobolev spaces, are defined by certain linear partial differential 
operators which may happen to depend on the PDEs under consider- 
ation. 

In the above example, for instance, the topology on the space X ob- 
viously depends on the specific PDE in (4.1). Thus the topological 
completion X in which the generalized solutions U are found accord- 
ing to (4.7), does again depend on the respective PDE. 



5. Pull-Back Partial Order 

Without loss of generality, [3] , we shall assume that all the posets con- 
sidered are without a minimum or maximum element. Various notions 
and results related to partial orders which are used in the sequel are 
presented in the Appendix. 

Given an equation (3.1), (3.2), we define on X the equivalence relation 
~t by 

(5.1) u^ T y T(u) = T(v) 

for u, v e X. In this way, by considering the quotient space 

(5.2) X T = X/^ T 
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we obtain the injective mapping 

(5.3) T^:X T — > Y 
denned by 

(5.4) X T 3Ui — ► T(u) G Y 

where u G U, that is, U is the ~t equivalence class of u in X T , while 
T(u) is defined by (3.2). 

At that stage, we can define a partial order <t on Xt as being the 
pull-back by the mapping T~ in (5.3) of the given partial order < on 
Y, namely 

(5.5) U < T V Tn(U) < T^(V) 

for U, V G X T . The effect of the above construction is that we obtain 
the order isomorphic embedding, or in short OIE 

(5.6) T~ : Xt — > Y 

As mentioned, without loss of generality we shall assume that the 
poset (AT*, <t) has no minimum or maximum. 

And now, we consider the order completions X* and Y* of Xt, and 
respectively, Y . 

For simplicity, we shall denote by < the partial orders both on A* 
and Y#. In fact, as seen in (A. 20), these partial orders are the usual 
inclusion relations C among subsets of A, respectively, of Y. 

Then according to Proposition A.l in the Appendix, we obtain the 
commutative diagram of OIE-s 
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Xn 



Y 



(5.7) 



X, 



c 



# 



c 



Y* 



Consequently, for U G X T and A G X*, we have in F # the relations 

(5.8) T*{<U\) = <T„(U)} 

(5.9) T#(A) = (T^A))* = sup y# { <T^(U)} \ UeA} 
6. Reformulation 

Now we can reformulate the problem of solving the general equations 
(3.1), (3.2) as follows. Given F G Y # , find necessary and sufficient 
conditions for the existence of A G X*, such that 

(6.1) T#(A) = F 



7. Solution 

We note that (5.7) gives the inclusions 

sup y# { T*(U) | U G X*, T*(U) C F } C 

C T#( sup x # { C/ | C/ G X#, T#(C/) C F } ) C 



(7.1) 



C T*{ inf x# { V | V G X#, F C T*(V) } ) C 



C inf y# { T#{V) \ V eX*, F C T # (V) } 



Indeed, the first and last inclusions follow from Lemma A.l in the 
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Appendix. As for the middle inclusion in (7.1), let U, V G X* be such 
that T*{U) C F C T*(V). Then T*{U) C T # (V), hence C/ < 
since T # is an OIE. It follows that 

sup x # { U | U G X*, T*{U) C F } < 

< inf x # { V | V G X#, F C T # (\/) } 

and the proof of (7.1) is completed. 

We note further that the above inequality T*{U) C F C T # (V) also 
implies 

sup r# { T # (f/) | [/ G X# T # (f/) C F } C F C 

(7.2) 

C inf y# { T # (V) | V G X#, F C T#(V) } 

Furthermore 

(7.3) {U eX* \ T#(U) C F } ^ 

since (A. 9) gives U = (p E X^, hence in view of (A. 4) - (A. 6) and 
(5.9) we have T*{U) = (T((f))) ul = <p ul = (<p u ) 1 = (Y*) 1 = <p C F. 

Returning now to the problem (6.1), we note that it is not trivial. 
Indeed, the OIE in (5.7), namely 

T* : X* — ► Y* 

need notbe surjective. Further T # need not preserve infima or suprema. 

However, in the next theorem we can obtain the following two general 
results : 

• a necessary and sufficient condition for the solvability of (6.1), 
and 

• the explicit expression of the solution, when it exists. 
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Theorem 7.1. 

Given F eY*. 

1) The equation 
(7.4) T*(A) = F 

has a solution A e X*, if and only if, see (7.1) 

sup y# { T*(U) \ U EX*, T*(U) CF} = 

(7.5) 

= inf y# { T#(V) \ V eX*, F C T*(V) } 

2) This solution is unique, whenever it exists, see (5.7). 

3) When it exists, the unique solution A e X* is given by 

A = sup x# { U e X* I T # ([7) C F } = 

(7.6) 

= inf x # {V eX* \ F C T#(V) } 

and, see (7.3) 

(7.7) { [/ e X# | T # (f/) C F }, {F6l*|F C T#(V) } ^ <f) 
Proof 

From (7.4) follows that AeX#, T#(A) C F, thus 

F = T # (A) C sup y# { T*{U) \ UeX#, T#{U) C F } 

Similarly we have 

inf y# { T # (V) | V G X#, F C T#{V) } C T # (A) = F 

thus (7.1) collapses to the seven equalities 

F = T#(A) = sup y# . . . = T#(sup x # . . . ) = 
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= T#(inf x # . . . ) = inf y# . . . = T*(A) = F 
Thus in particular we obtain (7.5). 

The injectivity of T* will give (7.6), while (7.7) follows from (7.3) and 
the fact that we can take V — A. 

Conversely, let us assume (7.5). Then (7.1) collapses to the three 
equalities 

sup y# ... = T # (sup x # . . . ) = T # (inf x # . . . ) = inf y # . . . 

thus in view of the corresponding collapsed version of (7.2), we can 
extend the above three equalities to the following four 

sup y# ... = T # (sup x # . . . ) = T # (inf x # . . . ) = inf y # ... = F 

And now the injectivity of T* will give (7.4) and (7.6), while (7.7) 
follows as above. 

□ 

The above existence result is of a "local" nature, since it refers to a 
solution of the equation (6.1) for one given right hand term F G F*. 
This result, however, can further be strengthened by the following 
"global" one which characterizes the solvability of (6.1) for all right 
hand terms F G Y*. Namely, we have, [3, pp. 190,191] 

Theorem 7.2. 

The following are equivalent 

(7.8) T*(X*) D Y 
and 

(7.9) T*{X*) = Y* 
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In each of these cases T# is an order isomorphism, or in short, 01, 
between X* and Y# . 

It is important to note the following four facts : 

• "Pull-back" type structures are customary when solving PDEs 
by functional analytic methods. Details in this regard are pre- 
sented in [3, chap. 12], while one well known classical example 
has been presented in section 4 above. 

• As shown in [3, chap. 13], and indicated in section 8 next, one 
can consider in (5.7) far more general partial orders than the 
"pull-back" type ones, and thus still obtain solutions by the 
order completion method for nonlinear systems of PDEs of type 

(1.1) , together with possibly associated initial and/or boundary 
value problems. 

• Nonlinear systems of PDEs with equations of type (1.1), together 
with possibly associated initial and/or boundary value problems 
can, in view of section 2 above, be dealt with as equations (3.1), 

(3.2) , and therefore (6.1). Consequently, the existence of solu- 
tions of such systems of PDEs follows from Theorems 7.1 and 
7.2. 

• The blanket, type independent, or universal regularity property 
of such solutions, namely, that they can be associated with usual 
measurable, or even Hausdorff continuous functions follows from 
arguments in [3, 1, 4-6]. 



8. Beyond "Pull-Back" Partial Orders 

As indicated here, and presented in full detail in [3, chap. 13], the 
order completion method in solving large classes of nonlinear systems 
of PDEs of the type (1.1) is not limited to the use of "pull-back" type 
partial orders in (45.7), (2.3) and (2.8). In fact, a large class of more 
general partial orders can be defined on the domains of the 
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respective PDEs, and stil obtain for them solutions in the correspond- 
ing order completions. 

As a brief indication in this regard, let us extend the approach in sec- 
tion 5 above as follows. Given any equation (3.1), (3.2) 

(8.1) T(A) = F 

where 



.2) T : X 



Y 



is any mapping, X is any nonvoid set, while (Y, <) is a poset, while 
F E Y is given, and A E X is the sought after solution. 

Further, let Z be any nonvoid set together with a surjective mapping, 
see (5.2) 

(8.3) A : Z — > X T 

Then we obtain the commutative diagram 



Y 



(8.4) 




Xn 



where T~ is given in (5.3), (5.4), while we have taken by definition 
(8.5) T\ = T^oA 

And now, we can define define the partial order <t, a on Z by 
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z — z', or 



Uz)<T x (z') 



(8.6) z < T , X z> 
In this case, we obtain 

(8.6) T~ is an OIE, T A is increasing, A is increasing and surjective 

and the procedures in sections 5-7 above can be applied to the map- 
ping 

(8.7) T\ : Z — ► Y 

instead of the initial mapping (3.2). 

Appendix 

We shortly present several notions and results used above. A related 
full presentation can be found in [3, Appendix, pp. 391-420]. 

Let (X, <) be a nonvoid poset without minimum or maximum. For 
a G X we denote 

(A.l) < a] = {x G X | x < a}, [a >= {x G X \ x > a} 

We define the mappings 

(A.2) X D A^A U = f] aeA [a > C X 

(A.3) X D A^A l = f] aeA < a] C X 

then for A C X we have 

(A.4) A u = X A 1 = X A = (f) 

(A. 5) A u = <fi A unbounded from above 
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(A. 6) A 1 = (f) •<=>- A unbounded from below 
Definition A.l. 

We call A C X a cut, if and only if 
(A.7) A ul = A 
and denote 

(A.8) X* = {A C X | A is a cut} C P(X) 

□ 

Clearly, (A.4) - (A. 6) imply 

(A.9) (f), x ex* 

therefore 

(A.10) X*^(f) 

Given A, B C X , we have 

(A.ll) A C B =>• A u D B u , A 1 D B l 

(A.12) A C A mZ , A C A iM 

(A. 13) = A u , A lul = A 1 

Consequently 
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V ACX : 

*) A ul e x* 

(A.14) * *) V B G X* : 

AC B =>■ A ul C £? 

therefore 

(A.15) X# = | A C X} 
Given x G X, we have 

(A. 16) {rr} u = [x >, {x} 1 =< x], [x > l =< x], < x] u = [x > 
(A.17) {x} ul =<x], {x} lu =[x> 
We denote for short 

Given A G X # , we have 

/ 3 a,beX : \ 

(A.18) 0^A^X^ 

\ < a] C A C < 6] / 
We shall use the embedding 

(A.19) X 3 x ^ x ul = x l =< x] G X # 

We define on X # the partial order 

(A.20) A < B ACB 
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Definition 2.1. 

Given two posets (X, <), (Y, <) and a mapping (p : X — > Y. We call 
ip an order isomorphic embedding, or in short, OIE, if and only if it is 
injective, and furthermore, for a, b G X we have 

a < b -<=/- </?(a) < 95(6) 

An OIE <£> is an order isomorphism, or in short, 01, if and only if it is 
bijective. 

□ 

The main result concerning order completion is given in, [2] : 

Theorem ( H M MacNeille, 1937 ) 

1) The poset (X*, <) is order complete. 

2) The embedding X X* in (A. 19) preserves infima and suprema, 
and it is an order isomorphic embedding, or OIE. 

3) For A E X*, we have the order density property of X in X*, 
namely 

A = sup x # {x l | x E X, x l C A} = 

(A.21) 

= inf x # {x l | x G X, AC x 1 } 

□ 

For A C X, we have 

(A.22) A ul = sup x# {x l | x G A} 

Given G AT#, with i G /, we have with the partial order in X# the 
relations 
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(A.23) sup, e/ Ai = inf {A e X# \ \J ieI A* C A} = (\J, l€l A t ) ul 

inf ieJ = sup {AeX# | A C f] ieI Ai} = (r\iei A i) ul = 

(A.24) 

= rw4 

Extending mappings to order completions 

Let (X, <), (Y, <) be two posets without minimum or maximum, and 
let 

(A.25) <p:X — > Y 

be any mapping. Our interest is to obtain an extension 

v # : x * — > r# 

For that, we first extend </? to a larger domain, as follows 
(A.26) y? # : PPO — ► ^ # 
where for A C X we define 

(A.27) y9#(A) = (¥>W = sup Y# {<^)]|x G yl} 

and for any mapping in (A.25), we obtain the commutative diagram 

X 3x i ► (p(x) G Y 

(A.28) 

V{X)3{x}\ k v #(x) = < V (x)]eY# 

Proposition A.l. 

1) The mapping </? # : V(X) — > Y* in (A. 36) is increasing, if on 
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V(X) we take the partial order defined by the usual inclusion "C". 

2) If the mapping tp : X — ► Y in (A. 35) is increasing, then the 
mapping : V(X) — > Y# in (A. 36) is an extension of it to X#, 
namely, we have the commutative diagram 

X 3 x i f. k ip(x) G Y 

(A29) 

X#3<x]\ ► Lp*(<x])= <Lp(x)]eY# 

3) If the mapping (p : X — > Y in (A. 25) is an OIE, then the mapping 
cp# : V(X) — ► Y# in (A.26) when restricted to X#, that is 

(A.30) p* :X* — > Y* 

as in (A.29), is also an OIE. 



Lemma A.l. 

Let in general /x : M — >■ N be an increasing mapping between two 
order complete posets, then for nonvoid E C M we have 

(A.31) /i(mi M E) < infjv /i(E) < sup N fi(E) < /x(sup M E) 

Proof 

Indeed, let a = inf M E e M. Then a < b, with b <E E. Hence 
/i(a) < /i(6), with b e E. Thus /i(a) < inf^v fi(E), and the first 
inequality is proved. 

The last inequality is obtained in a similar manner, while the middle 
inequality is trivial. 

□ 
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